Naval  Ocean  Research  and 
Development  Activity 

NSTl.  Mississippi  39529 


$8 


!»■ 


'JK 


:  Vf 


■A  ;m  m 


NORDA  Technical  Note  272 


CHORDS:  A  New  Temperature  or 
Sound  Speed  Profile  Thinning  Algorithm 


SECURITY  CLASSIFICATION  OF  THU  PAPE  (Whm*  Dolt  Knitted) 

REPORT  DOCUMENTATION  PAGE 


|T  PERFORMING  ORGANIZATION  NAME  AND  ADDRESS 

Naval  Ocean  Research  &  Development  Activity 
Ocean  Acoustics  and  Technology  Directorate 
NSTL.  Mississippi  39529 


[T  REPORT  NUMBER 

NORDA  Tecnnical  Note  272 


4.  TITLE  fend  Submit) 

CHORDS:  A  New  Temperature  or  Sound  Speed 
Profile  Thinning  Algorithm 


”  It.  OOVT  ACCESSION  n£  t  RECIFE* 

Ab 

I.  TYRE  OF 


7.  AUTHONf#; 

George  Kerr 


II.  CONTROLLING  OFFICE  NAME  AND  ADDRESS 

Same 


14.  MONITORING  AOENltv  NAME  A  ADORESSfl/  dlllettnt  Tiorn  Conlrolllnt  Olll cb) 


14.  DISTRIBUTION  STATEMENT  (el  (hit  Iftporl) 

Approved  for  Public  Release 
Distribution  Unlimited 


READ  ZMITRUCTtOMB 
BEFORE  COttPLETWO  FORM 

FE.CI  PI  EN  T'S  CATALOG  NUMBER 


I.  TYPB  OF  REPORT  *  PERIOD  COVERED 

Final 


i.  pi 


RIM 


•garrwff  grwTHriroaiw?g 


10.  PRO 


wimMmw™ 

P.E.  940101 


It.  REPORT  DATE 


II.  NUMBER  OF  PAOBI 


30 


II.  ISCURITY  CL  All.  (el  mil  1*553“™“ 

UNCLASSIFIED 

TTT '  c  AfidNTa&ffNA'rATngy  • 


17.  DISTRIBUTION  IT  ATEMI1  NT  (a I  Iht  tbtlrtcl  intend  In  tlock  SO,  II  dlllttml  tnm  Ktpotl) 


t|.  SUPPLEMENTARY  NOTES 


It.  KEY  WORDS  rConl In 


lory  aiw  Ittnltty  by  block  number) 


profile  thinning 
sound  speed 


temperature 
digital  data 
profile  filtering 


20 


ABSTRACT  (Continue  on  rovoroo  aldo  H  nocottory  md  Identity  fry  k/oek  numhof) 

.  A  new  sound  velocity  and  temperature  profile  thinning  algorithm  especially 
designed  for  use  with  high  vertical  resolution  (1  meter)  conductivity- 
temperature-depth  (CTD)  or  expendable  bathythermograph  (XBT)  digital  data 
is  described.  The  thinning  operation  can  be  halted  once  the  desired  number 
of  output  points  is  obtained  or  alternatively  when  a  user-supplied  tolerance 
value  is  reached.  The  new  algorithm  has  the  advantage  of  scanning  the  entire 
profile  during  each  iteration.  As  a  result,  unlike  other  thinning  algorithms 
examined,  profile  curvature  at  depth  is  retained.  _ 


DD  lETn  1473 


EDITION  OF  I  NOV  *1  II  OBSOLETE 

S/N  0103- LF.  014-  6401 


SECURITY  CLASH  Fie  A1 


ILLMHEL 

iification  or  this 


page  dhisN  Ms  Msisrssj 


ABSTRACT 


- >  A  new  sound  speed  and  temperature  profile  thinning  algorithm  especially 

designed  for  use  with  high  vertical  resolution  (1  meter)  conductivity -temperature- 
depth  (CTD)  or  expendable  bathythermograph  (XBT)  digital  data  Is  described.  The 
thinning  operation  can  be  halted  once  the  desired  number  of  output  points  Is 
obtained  or  alternatively  when  a  user-supplied  tolerance  value  Is  reached.  The  new 
algorithm  has  the  advantage  of  scanning  the  entire  profile  during  each  Iteration.  As 
a  result,  unlike  other  thinning  algorithms  examined,  profile  curvature  at  depth  Is 
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CHORDS:  A  MEW  TEMPERATURE  OR  SOUND  SPEED  PROFILE  THINNING  ALGORITHM 


I.  INTRODUCTION 

Digitally  recording  oceanographic  sensors  have  the  capability  of  gathering 
hundreds  or  even  thousands  of  data  points  In  a  single  cast.  In  contrast,  present 
acoustic  models  are  capable  of  utilizing  only  a  small  number  (on  the  order  of  30)  of 
these  points.  It  Is,  therefore.  Important  to  have  a  reliable  method  of  thinning  the 
digital  data  to  the  number  of  points  acceptable  to  the  acoustic  models.  In  addition 
to  speed,  automated  methods  offer  the  advantage  of  uniformity  over  "by  hand" 
operations. 

Two  automated  thinning  algorithms  (NORDA's  FRITZ  and  an  algorithm  obtained  from 
the  SACLANT  Center)  were  examined  to  determine  their  performance  with  T7-XBT  profiles 
with  1  meter  (m)  resolution  from  the  surface  to  100  m  and  10  m  resolution  from  100  m 
to  the  bottom  of  the  XBT  (approximately  800  m).  Both  algorithms  were  similar  In 
function  with  the  SACLANT  Center  algorithm  being  somewhat  more  sophisticated.  These 
algorithms  contained  a  particular  flaw  that  made  them  unacceptable.  Both  algorithms 
started  processing  at  at  the  surface  and  worked  down  the  profile.  As  a  result  the 
bottom  portion  of  the  temperature  profile  was  often  not  reached.  In  other  words,  the 
desired  number  of  output  points  were  found  before  the  entire  profile  was  operated  on 
by  the  algorithms.  As  a  result,  the  curvature  In  the  bottom  portion  of  the  measured 
profile  was  removed  from  the  thinned  profile. 

Figure  1  Illustrates  a  complex  staircased  sound  speed  profile  of  approximately 
1  m  resolution.  Figure  2  Illustrates  the  same  sound  speed  profile  thinned  to  30 
points  using  the  SMOOTH  thinning  algorithm  (similar  In  methodology  to  the  algorithms 
mentioned  above)  Incorporated  Into  PRISM,  an  acoustic  analysis  model.  Notice  the 
change  In  the  sound  speed  gradient  In  the  layer,  the  loss  of  some  detail  In  the 
staircase  near  100  m,  and  the  loss  of  profile  curvature  below  400  m.  The  loss  of 
curvature  In  this  case  resulted  In  the  movement  of  the  deep  sound  channel  axis  from 
Its  actual  depth  of  near  800  m  to  a  depth  of  approximately  1000  m. 

Because  of  this  type  of  performance,  It  was  decided  to  design  a  thinning 
algorithm  that  would  Insure  that  the  entire  profile  was  considered  during  each  point 
selection  Iteration  of  the  algorithm.  The  result  Is  a  FORTRAN  subroutine  named 
CHORDS. 

II.  PROGRAM  GENERAL  DESCRIPTION 

As  the  name  Implies,  CHORDS  fits  straight  line  segments  (chords)  to  temperature 
or  sound  speed  profiles  (a  version  of  CHORDS  Is  available  for  both).  The  chords  are 
selected  so  that  the  maximum  difference  between  any  new  chord  and  the  profile  seg¬ 
ment  It  encompasses  Is  diminished  as  each  chord  Is  added  to  the  chord  set  defining 
the  entire  profile. 

The  CHORD  process  begins  with  the  selection  of  an  Initial  set  of  significant 
points.  For  a  temperature  profile  the  Initial  significant  points  are  the  first 
(surface.)  and  last  (bottom)  points  of  the  profile,  and  a  point  selected  to  represent 
the  mixed  layer  depth.  The  mixed  layc.-  depth  point  Is  selected  as  the  (n-l)th  point 
where  the  nth  point  temperature  differs  by  at  least  ±.2°C  from  the  surface  tempera¬ 
ture.  The  three  resulting  significant  points  then  define  an  Initial  set  of  two 
chords . 
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For  a  sound  speed  profile,  the  significant  points  are  the  first  and  last  pro¬ 
file  points  and  points  defining  the  layer  depth  and  depth  of  the  deep  sound  channel 
axis.  The  deep  sound  channel  axis  depth  is  defined  as  the  depth  at  which  the  abso¬ 
lute  sound  speed  mlnlmun  occurs.  The  layer  depth  Is  defined  as  the  depth  above  the 
deep  sound  channel  axis  depth  where  the  local  maximum  sound  speed  occurs.  Thus  for 
sound  speed  profiles  an  Initial  set  of  up  to  three  chords  Is  possible.  Figure  3 
Illustrates  an  Initial  point  and  chord  set. 

The  algorithm  then  treats  each  chord  and  associated  profile  segment  In  the 
following  manner: 

A— The  chord  and  corresponding  profile  segment  are  rotated  and  translated 
(Figure  4)  so  that  the  chord  becomes  the  new  temperature  (sound  speed)  axis. 

B— The  absolute  maximum  of  the  transformed  profile  segment  point  set  Is  then 
located  (Figure  5).  This  point,  defining  the  maximum  absolute  difference  between  the 
profile  segment  and  Its  associated  chord,  becomes  a  new  significant  point  candidate. 

After  all  chords  are  processed,  the  candidate  significant  points  are  compared. 
The  candidate  point  having  the  largest  absolute  difference  between  Itself  and  Its 
associated  chord  of  all  such  pairs  Is  selected  as  a  new  significant  point.  The  new 
point  divides  an  existing  profile  segment  Into  two  parts  (Figure  5),  thereby  creat¬ 
ing  two  chords  In  place  of  a  single  previously  existing  chord.  By  retaining  the 
chord  to  profile  segment  maximum  absolute  differences  of  all  chord-profile  segment 
pairs,  further  Iterations  will  only  require  two  chord-profile  segment  transforma¬ 
tions  and  maximum  difference  calculations  (one  Iteration  for  each  of  the  two  newly 
created  chords). 

The  above  chord-splitting  process  will  continue  until  the  number  of  desired 
points  Is  found,  or  until  the  maximum  allowable  absolute  difference  between  any 
chord  and  Its  encompassed  profile  segment  Is  less  than  a  user-defined  tolerance 
limit.  In  tests  of  the  algorithm's  performance,  a  metric  unit  tolerance  value  of 
0.05  for  temperature  and  0.1  for  sound  speed  were  found  to  adequately  reproduce 
profile  shapes  and  yet  reduce  the  number  of  profile  points  by  approximately  a  factor 
of  10. 

III.  PROGRAM  EXAMPLE 

Figure  6  shows  the  test  sound  speed  profile  (Figure  1),  addressed  previously, 
thinned  to  30  points  using  the  CHORDS  algorithm.  Notice  how  the  In-layer  gradient  of 
the  original  profile  has  been  retained  along  with  much  of  the  staircase  structure. 
Notice  also  how  the  curvature  of  the  profile  below  400  m  has  been  retained.  Here  the 
shape  of  the  original  profile  Is  well -matched  In  the  vicinity  of  the  deep  sound 
channel ,  and  the  depth  of  the  deep  sound  channel  axis  has  been  retained. 

Appendix  A  gives  a  more  detailed  description  of  the  subroutine  In  addition  to 
Instructions  for  Its  use.  Appendix  B  contains  a  listing  of  the  program.  Appendix  C 
contains  a  description  of  the  subroutines,  functions,  and  variable  names  used  In  the 
program. 

IV.  CONCLUSIONS 

CHORDS  has  been  found  to  produce  a  more  accurate  representation  of  both 
high-resolution  temperature  and  sound  speed  profiles  than  several  other  profile 
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Figure  4.  Translation  and  rotation  of  a  chord  and  corresponding  profile  segment 
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thinning  algorithms  In  a  limited  r.irnbe*  of  test  cases.  Since  profiles  covering  the 
extremes  from  fairly  smooth  to  highly  variable  (the  staircased  profile  presented  in 
this  report)  were  tested,  it  is  reasonable  to  assume  that  CHORDS  will  also  perform 
equally  as  well  for  profiles  between  these-  extremes.  For  low-frequency  applications 
It  Is  possible  that  other  algorithms  may  be  modified  to  enable  them  to  better  repre¬ 
sent  the  deeper  portions  of  the  profiles  being  thinned.  However,  for  high-frequency 
work  or  high-resolution  highly  variable  input  profiles,  simple  fixes  may  not  be 
possible. 

It  Is  suggested  that,  potential  users  of  CHORDS  pa}  particular  attention  to 
their  own  applications  and  test  CHORDS  accordingly.  The  tolerance  values  suggested 
and  the  number  of  output  points  required  may  not  be  optimal  for  all  applications. 

V.  RECOMMENDATIONS 

It  Is  recommended  that  CHORDS  be  seriously  considered  In  all  applications 
(i.e.,  NORDAPS,  TOPS  Input  to  SHARPS)  where  a  reduction  in  the  number  of  sound  speed 
or  temperature  profile  points  is  required. 

The  CHORDS  concept  may  be  easily  applied  to  the  thinning  problems  (i.e., 
density  profiles,  salinity  profiles,  high-resolution  numerical  model  output,  etc.). 
These  other  applications  should  be  given  additional  thought  and  testing. 

With  the  capability  now  available  to  capture  and  store  high  resolution  (I.e.,  1 
meter  temperature  and  sound  speed)  Information,  It  Is  Important  to  know  the  effect 
of  thinning  this  data  to  some  arbitrary  lower  resolution.  It  Is  therefore 
recommended  that  studies  be  made  In  areas  where  results  of  such  thinning  may  have 
Important  consequences.  One  such  area  for  study  would  be  the  relationship  between 
transmission  loss  and  the  number  of  retained  sound  speed  profile  points.  It  Is 
expected  that  such  a  study  would  result  In  the  definition  of  the  optimal  number  of 
profile  points  required  to  obtain  accurate  transmission  loss  results  In  various 
ocean  areas.  In  addition  to  being  area  dependent,  results  are  expected  to  depend  on 
the  closeness  of  the  profile  match  at  various  profile  depths  as  related  to  the 
effected  modes  of  sound  propagation. 


APPENDIX  A:  CHORDS  DETAILED  DESCRIPTION 


I.  INPUT  AND  OUTPUT 

CHORDS  Is  a  FORTRAN  V  callable  subroutine  written  on  a  Digital  Equipment 
Corporation  VAX  11/780  minicomputer.  The  subroutine  reduces  the  number  of  Input 
sound  speed  profile  points  using  the  desired  number  of  output  points  or  a  tolerance 
limit  to  cease  processing.  The  tolerance  limit  refers  to  the  maximum  acceptable 
difference  between  a  chord  and  Its  corresponding  profile  segment.  By  passing  a 
positive  tolerance  limit  to  the  subroutine,  output  profiles  of  varying  numbers  of 
points  will  result.  The  number  of  output  points  will  depend  on  the  tolerance  limit 
selected  and  the  complexity  of  the  profile  being  thinned.  Tolerance  limit  values  of 
0.1  for  sound  speed  profiles  and  0.05  for  temperature  profiles  have  reproduced  test 
profile  shapes  sufficient  for  low  frequency  acoustic  applications  while  reducing  the 
number  of  profile  points  by  a  factor  of  10.  A  zero  value  of  the  tolerance  limit  will 
force  the  subroutine  to  halt  the  thinning  operation  when  the  total  number  of  desired 
output  points  has  been  reached. 

There  are  three  arguments  In  the  subroutine  call: 

CALL  CHORDS  (NIN  ,  N0UT,  T0L) ,  where 

NIN  *  the  user-supplied  number  of  Input  profile  points,  N0UT  »  the  number  of 
output  profile  points  determined  by  the  algorithm,  and 

T0L  *  the  user-supplied  tolerance  value. 

The  Input  profile  is  passed  to  the  subroutine  through  the  following  named 
common  area: 

COMMON  /DATA/  D(ASIZ) ,S( ASIZ) ,  where 

D  Is  the  Input  depth  array,  and 

S  Is  the  Input  sound  speed  array. 

The  output  profile  Is  passed  through  the  following  named  common  area: 

COMMON  /DATAO/  OR(MNPTS) ,SR(MNPTS) ,  where 

DR  Is  the  output  depth  array,  and 

SR  is  the  output  sound  speed  array. 

The  program  Is  set  up  to  handle  profiles  consisting  of  a  maximum  of  1000 
points.  This  number  can  be  changed  by  modifying  the  parameter  ASIZ,  In  all 
subroutines  In  which  this  parameter  appears,  to  equal  the  maximum  number  of  Input 
profile  points  expected. 

The  maximum  number  of  output  profile  points  is  presently  set  to  30.  This  number 
may  be  changed  by  modifying  the  parameter  MNPTS  appearing  In  the  parameter  statement 
In  subroutine  CHORDS.  If  enough  profile  points  are  not  supplied  to  the  subroutine,  a 
zero  value  will  be  returned  for  the  number  of  output  points  (NR). 
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II.  DETAILED  FUNCTIONAL  DESCRIPTION 


In  the  following  discussion  the  numbers  In  parentheses  refer  to  the  program 
line  numbers  Involved.  A  program  listing  can  be  found  In  Appendix  B.  A  complete 
variable  list  and  externals  called  can  be  found  In  Appendix  C.  A  flow  chart 
Indicating  the  relationship  between  the  subroutines  and  functions  called  can  be 
found  In  Figure  Al.  The  subroutine  discussed  Is  designed  for  use  with  sound  speed 
profiles.  The  version  designed  for  temperature  profiles  differs  only  In  the 
selection  of  the  Initial  set  of  chord  points. 

As  outlined  In  Figure  A'i,  CHORDS  first  chooses  an  Initial  set  of  chords.  For 
sound  speed  profiles  the  Initial  set  may  be  composed  of  from  two  to  four  points  or 
from  one  to  three  chords. 

The  Initial  set  of  points  will  always  contain  the  surface  (depth  *  0)  point 
(93)  and  the  last  profile  point  (981.  If  the  first  profile  point  Is  not  at  the 
surface,  then  a  surface  value  of  s  d  speed  Is  determined  by  a  call  to  subroutine 
INTERP  (67)  which  uses  the  first  twc  profile  points  to  extrapolate  sound  speed  to 
the  surface.  The  first  profile  point  Is  replaced  by  the  surface  point.  Other  Initial 
points  could  be  the  point  defining  the  layer  depth  (84-88)  and  the  point  defining 
the  deep  sound  channel  axis  depth  (75-79).  The  deep  sound  channel  axis  depth  Is 
selected  as  that  depth  where  the  minimum  sound  speed  Is  found.  The  layer  depth  Is 
defined  as  that  depth  where  the  above  axis  maximum  sound  speed  Is  found.  These 
additional  Initial  points  will  not  be  used  If  they  coincide  with  the  surface  or  last 
profile  points. 

The  subroutine  designed  for  use  with  temperature  profiles  differs  from  the 
above  In  that  the  only  addition  Initial  point  which  might  be  added  to  the  surface 
and  bottom  points  Is  a  point  defining  the  mixed  layer  depth.  This  depth  Is  defined 
as  that  depth  below  which  the  temperature  differs  from  the  surface  temperature  by 
±0.2°C. 

For  each  Initial  chord,  the  encompassed  profile  segment  point  farthest  away 
from  the  chord  Is  located  through  a  call  to  subroutine  CHORD  (144).  Subroutine  CHORD 
calls  subroutine  RTRAN  which  translates  and  rotates  the  profile  segment  points  so 
that  the  chord  becomes  the  new  X  axis.  Subroutine  RTRAN  uses  function  ALPHA  to 
determine  the  angle  of  rotation.  The  profile  segment  point  lying  farthest  from  the 
now  horizontal  chord  can  be  easily  found  by  locating  the  point  of  largest  absolute 
magnitude.  This  magnitude  and  the  corresponding  profile  point  number  are  returned 
and  stored  in  arrays  DIFF  and  NUM  respectively. 

The  resulting  difference  set  (array  DIFF)  Is  then  examined  to  determine  which 
element  possesses  the  largest  chord  to  profile  segment  difference  (188-191).  The 
profile  point  at  which  this  maximum  difference  Is  found  (stored  In  array  NUM) 
becomes  an  end  point  for  two  new  chords  (235-239).  For  one  new  chord,  the  point  will 
be  the  first  chord  point.  For  the  other  new  chord,  the  point  will  be  the  last  chord 
point.  The  chord  counter  LINCNT  Is  Incremented  by  one  (196)  since  the  two  new  chords 
replace  an  existing  chord,  thereby  Increasing  the  total  number  of  chords  by  one. 

The  final  step  for  this  first  new  chord,  determined  from  the  Initial  chord  set, 
Is  to  store  the  new  chord  end  points  (array  PTPAIR) ,  and  make  room  In  the  difference 
set  array  (DIFF)  and  the  point  of  maximum  difference  array  (NUM)  for  the  yet  to  be 
determined  values  of  these  variables  for  the  two  new  chords.  The  section  of  code 
which  handles  the  array  Insertion  depends  on  where  In  the  chord  sequence  the  new 
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chord  Is  found  (217-226).  The  section  beginning  at  program  line  227  and  ending  at 
program  line  251  deals  with  the  Insertion  of  the  two  new  chords  In  the  arrays  when 
the  chord  number  counter  J  and  the  new  chord  sequence  number  (MAXLIN)  match  (221). 
The  variable  MAXLIN  Is  used  to  determine  In  which  original  chord  the  new  two  chords 
are  located.  Variable  LMAX  Is  equated  to  MAXLIN  and  variable  MAX  Is  made  equal  to 
MAXLIN  plus  one  to  retain  the  new  chord  sequence  numbers  for  the  next  Iteration.  The 
program  section  beginning  at  line  256  and  extending  to  line  268  bumps  the  existing 
arrays  down  after  Insertion  of  the  new  data.  The  program  section  beginning  at  line 
272  and  ending  at  line  280  Inserts  the  last  or  bottom  chord  Information  In  the 
arrays. 

In  the  following  Iterations,  the  only  chord  to  profile  segment  maximum  differ¬ 
ences  that  need  to  be  calculated  are  for  the  two  chords  defined  In  the  previous 
Interatlon.  These  differences  are  found  through  the  call  to  CHORD  found  at  program 
line  177.  The  CHORD  call  Is  triggered  by  a  match  of  the  chord  number  counter  J  to 
either  variable  LMAX  or  MAX  (160).  In  the  loop  beginning  at  program  line  149  the 
maximum  chord  to  profile  segment  difference  Is  again  located  for  the  chord  set  now 
Including  the  two  new  chords  (188-191).  The  chord-splitting  process  Is  repeated 
until  either  the  total  number  of  desired  points  have  been  located  (206)  or  until  the 
tolerance  limit  Is  reached  (200). 

Processing  stops  when  either  of  the  two  above  conditions  are  met.  The  resulting 
profile  (arrays  DR  and  SR),  except  for  the  last  profile  point,  Is  generated  from  the 
first  elements  of  the  chord  end  point  array  (PTPAIR).  The  last  profile  point  Is 
defined  by  the  last  of  the  second  elements  of  the  same  array. 


Figure  Al.  Functional  relationship  between  the  main  calling  program  and  CHORDS 
subroutines  and  functions 
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Figure  A2.  Logical  flow  diagram  for  the  CHORDS  point  selection  algorithm 
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APPENDIX  0:  SUBROUTINE  LISTING 


0001 

0002  C 

0003  C 

0004  C 

0005  C 

0006  C 

0007  C 

0008  C 

0009  C 

0010  C 

0011  C 

0012  C 

0013  C 

0014  C 

0015  C 

0016  C 

0017  C 

0018  C 

0019  C 

0020  C 

0021  C 

0022  C 

0023  C 

0024  C 

0025  C 

0026  C 

0027  C 

0028  C 

0029  C 

0030  C 

0031  C 

0032  C 

0033  C 

0034  C 

0035  C 

0036 

0037  C 

0038 

0039 

0040 

0041 

0042 

0043 

0044 

0045  ' 

0046  C 

0047  C 

0048  C 

0049 


SUBROUTINE  CHORDS! N IN, NR ,CTDLIM) 

THIS  SUBROUTINE  PERFORMS  A  THINNING  OF  SMOOTHED 
OR  EDITED  SOUND  SPEED  PROFILES.  FOR  A  FILTERING 
ALGORITHM  USED  FOR  CTD  1  METER  DATA  CONTACT  GEORGE  KERR 
NORDA  CODE  323,  NSTL,  MS  39529,  OR  EUGENE 
M0L1NELLI,  PSI ,  MCLEAN,  VA.( AUTHOR  OF  THE  INITIAL 
FILTERING  ALGORITHM). 

THIS  SUBROUTINE  IS  CURRENTLY  SET  UP  TO  THIN  SOUND  SPEED 
PROFILES  CONSISTING  OF  A  MAXIMUM  OF  1000  POINTS  TO  A 
PROFILE  CONSISTING  OF  EXACTLY  30  POINTS.  THE  MINIMUM 
NUMBER  OF  INPUT  PROFILE  POINTS  IS  32. 

TO  CHANGE  THESE  LIMITS  PARAMETERS  MPNTS(DEFINING  THE  NUMBER 
OF  POINTS  TO  FILTER  TO),  AND  ASIZ(DEFINING  THE  MAXIMUM  NUMBER 
OF  INPUT  PROFILE  POINTS),  MUST  BE  CHANGED  IN  ALL  SUBROUTINES 
IN  WHICH  THESE  PARAMETERS  APPEAR. 


THE  MAIN  PROGRAM  CALLING  THIS  ROUTINE  MUST  CONTAIN 
THE  COMMON  AREAS  /DATA/,  AND  /DAT AO/  TO  HOLD  THE  RAW  INPUT 
DEPTHS  (D)  AND  SOUND  SPEEDS  (S),  AND  THE  OUTPUT  DEPTHS  (DR) 
AND 

SOUND  SPEEDS  (SR)  RESPECTIVELY. 

THE  ONLY  ARGUMENTS  PASSED  TO  THE  SUBROUTINE  ARE  THE  NUMBER 
OF  INPUT  DATA  POINTS  (NIN)  AND  THE  DESIRED  TOLERANCE  LIMIT 
(CTDLIM) .  NR,  RETURNED  BY  THE  SUBROUTINE  AS  A  CALLING 
ARGUMENT, 

WAS  INCLUDED  TO  RETURN  THE  NUMBER  OF  POINTS  THINNED  TO 
IF  THE  TOLERANCE  LIMIT  RATHER  THAN  THE  NUMBER  OF  POINTS  TO 
THIN 

TO  WAS  EMPLOYED  TO  HALT  THE  THINNING  OPERATION. 

SUBROUTINE  CHORDS  DEVISED  AND  PROGRAMMED  BY  GEORGE  KERR 
NORDA  CODE  223,  NSTL,  MS  39529. ( 2MAY83) 

PARAMETER  ASIZ-1000.MPNTS-30 

COMMON  /DATA/  D( ASIZ) ,S( ASIZ) 

COMMON  /DATAO/  DR(MPNTS) ,SR(MPNTS) 

COMMON  /IPOINT/  IBEG.IEND 

INTEGER  PTPAIR(MPNTS,2), TEMP  1 , TEMP 2 , TEMP 3 , T EMP 4 

REAL  INTERP 

DIMENSION  DIFF(MPNTS) .NUM(MPNTS) 

DATA  DIFFLT/-999S9./ 

NOLIN-MPNTS-1 

ZERO  ARRAYS 

DO  9999  U-l.MPNTS 


pf»«VIOU»  PAO* 
It  BLANK 
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0050 

0051 

0052 

0053 

0054 

0055 

0056 

0057 

0058 

0059 

0060 

0061 

0062 

0063 

0064 

0065 

0066 

0067 

0068 

0069 

0070 

0071 

0072 

0073 

0074 

0075 

0076 

0077 

0078 

0079 

0080 

0081 

0082 

0083 

0084 

0085 

0086 

0087 

0088 

0089 

0090 

0091 

0092 

0093 

0094 

0095 

0096 

0097 

0098 

0099 

0100 

0101 

0102 

0103 


DIFF( J)“0. 

NUM(J )=0 
9999  CONTINUE 
NR-0 
C 

C  DEFINE  THE  MINIMUM  NUMBER  OF  POINTS  ALLOWED 

C 

MINPTS-NOLIN+3 

C 

C  CHECK  TO  SEE  IF  THE  PROFILE  HAS  SUFFICIENT  POINTS 

C 

IF(NIN.LT.MINPTS)  RETURN 
C 

C  IF  THE  PROFILE  DOES  NOT  HAVE  A  SURFACE  VALUE 

C  EXTRAPOLATE  TO  FIND  A  SURFACE  VALUE 

C 

IF(D( 1)  .EQ.O.Q)  GO  TO  2500 
S(1)«INTERP(D(1),S(1),D(2),$(2),0.0) 

2500  D(1)«0.0 
SSMIN-99999. 

SSMAX--99999. 

C 

C  FIND  THE  DEEP  SOUND  CHANNEL  AXIS  (DEFINED  HERE  AS 

C  THE  ABSOLUTE  MINIMUM  OF  ALL  SOUND  SPEEDS  INPUT) 

C 

DO  3000  J-l.NIN 

IF(S( J) . GE.SSMIN)  GO  TO  3000 

MINPT-J 

SSMIN-S(J) 

3000  CONTINUE 
C 

C  FIND  THE  LAYER  DEPTH  (DEFINED  HERE  AS  THE  ABSOLUTE 

C  SOUND  SPEED  MAXIMUM  ABOVE  THE  DEEP  SOUND  CHANNEL  AXIS  DEPTH) 

C 

DO  3001  J-MINPT.1,-1 
IF(S(J) .LE.SSMAX)  GO  TO  3001 
MAXPT-J 
SSMAX-S(J) 

3001  CONTINUE 
C 

C  RETAIN  THE  SURFACE  POINT  AS 

C  THE  FIRST  POINT  ON  THE  FIRST  CHORD 

C 

PTPAIR(1,1)«1 

C 

C  RETAIN  THE  BOTTOM  POINT  AS 

C  THE  LAST  POINT  ON  THE  FIRST  CHORD 

C 

PTPAIR(1,2)-NIN 

C 

C  INITIALIZE  THE  NUMBER  OF  CHORDS  FOUND 

C 

LINCNT-1 

C 
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0104 

0105 

0106 

0107 

0108 

0109 

0110 

0111 

0112 

0113 

0114 

0115 

0116 

0117 

0118 

0119 

0120 

0121 

0122 

0123 

0124 

0125 

0126 

0127 

0128 

0129 

0130 

0131 

0132 

0133 

0134 

0135 

0136 

0137 

0138 

0139 

0140 

0141 

0142 

0143 

0144 

0145 

0146 

0147 

0148 

0149 

0150 

0151 

0152 

0153 

0154 

0155 

0156 

0157 


C  CHECK  TO  INSURE  THAT  THE  LAYER  DEPTH  POINT  IS 

C  NOT  THE  SAME  AS  THE  FIRST  OR  LAST  PROFILE  POINT 

C 

lF(MAXPT.EQ.l.QR.MAXPT.EQ.NIN)  GO  TO  83 
C 

C  INSERT  THE  LAYER  DEPTH  DEFINED  CHORDS  IN  THE  CHORD 

C  END  POINT  ARRAY 

C 

PTPA1R(LINCNT,2)«MAXPT 
LINCNT-L1NCNT+1 
PTPAIR(LINCNT,1)-MAXPT 
PTPAIR(LINCNT,2)*NIN 
83  CONTINUE 
C 

C  MAKE  SURE  THE  DEEP  SOUND  CHANNEL  AXIS  DEPTH 

C  POINT  IS  NOT  THE  SAME  AS  THE  FIRST,  LAST,  OR 

C  LAYER  DEPTH  POINT 

C 

IF(MINPT.EQ.l.QR.MINPT.EQ.NIN.OR.MINPT.EQ.MAXPT)  GO  TO  80 
C 

C  INSERT  THE  DEEP  SOUND  CHANNEL  AXIS  DEFINED  CHORDS 

C  IN  THE  CHORD  END  POINT  ARRAY 

C 

PTPAIR(LINCNT,2)-MINPT 
LINCNT-LINCNT+1 
PTPAIR(LINCNT,1)«MINPT 
PTPAIR(LINCNT,2)-NIN 
80  CONTINUE 

DO  85  J-l.LINCNT 
C 

C  CHOOSE  CHORD  END  POINTS  SO  THAT  ALL  POINTS 

C  ALONG  THE  CHORD  CAN  BE  IDENTIFIED. 

C 

IBEG-PTPAIR(J.l) 

IEND-PTPAIR(U,2) 

C 

C  FOR  EACH  CHORD  DETERMINE  THE  MAXIMUM  DEVIATION  FRuM 

C  THE  ROTATED  AXIS  AND  THE  NUMBER  OF  THE  POINT  AT  WHICH 

C  THIS  OCCURS. 

C 

CALL  CHORD(DIFF( J ) ,NUM(U)) 

85  CONTINUE 
NMAX-0 
LMAX-0 

110  DIFMAX-99999. 

DO  200  J-l.LINCNT 
C 

C  LOOK  FOR  NEWLY  DEFINED  CHORDS  FOR  WHICH  MAXIMUM 

C  DIFFERENCES  HAVE  NOT  BEEN  FOUND.  NEW  CHORDS  WILL 

•C  BE  DEFINED  BY  A  PREVIOUS  COMPARISON  OF  ALL  CHORD 

C  DIFFERENCES.  THE  CURVE  HAVING  THE  LARGEST  DIFFERENCE 

C  FROM  ITS  CHORD  WILL  DEFINE  2  NEW  CHORDS.  ONLY 

C  THESE  TWO  NEW  CHORDS  NEED  TO  HAVE  DIFFERENCES 

C  CALCULATED.  DIFFERENCES  FOR  THE  REMAINING  CHORDS 
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0158 

C 

HAVE  ALREADY  BEEN  FOUND. 

0159 

C 

0160 

IF(J.EQ.NMAX.OR.J .EQ.LMAX)  GO  TO  154 

0161 

GO  TO  155 

0162 

C 

0163 

c 

CALCULATE  NEW  CHORD  DIFFERENCES. 

0164 

c 

0165 

154 

IBEG«PTPAIR(J,1) 

0166 

IEND«PTPAIR(J,2) 

0167 

C 

0168 

C 

IF  THE  TWO  END  POINTS  OF  THE  CHORD  ARE  ADJACENT 

0169 

C 

NO  FURTHER  DIVISION  OF  THIS  CHORD  IS  POSSIBLE. 

0170 

C 

THEREFORE,  DO  NOT  CALCULATE  A  DIFFERENCE. 

0171 

C 

0172 

IF((IEND-IBEG).EQ.l)  GO  TO  160 

0173 

c 

0174 

c 

CALCULATE  THE  DIFFERENCE  AND  FIND  THE  POINT  NUMBER  FOR 

0175 

c 

THE  NEWLY  FOUND  CHORD. 

0176 

c 

0177 

CALL  CHORD (D IFF (J).NUM(J)) 

0178 

GO  TO  155 

0179 

c 

0180 

c 

FOR  ADJACENT  CHORD  END  POINTS 

0181 

c 

SET  THIS  CHORD  DIFFERENCE  EQUAL  TO  THE  DEFAULT  VALUE 

0182 

c 

0183 

160 

DIFF(J)«CTDLIM 

0184 

C 

0185 

c 

LOOK  FOR  THE  MAXIMUM  DIFFERENCE  AMONG  ALL 

0186 

c 

THE  CHORD  DIFFERENCES. 

0187 

c 

0188 

155 

IF ( D IFF < J ) . LT .DIFMAX)  GO  TO  200 

0189 

DIFMAX-DIFF(J) 

0190 

MAXLIN-J 

0191 

MAXPT-NUM(J) 

0192 

200 

CONTINUE 

0193 

c 

0194 

c 

INCREASE  CHORD  COUNTER 

0195 

c 

0196 

LINCNT-LINCNT+1 

0197 

c 

0198 

c 

HAS  THE  TOLERANCE  LIMIT  BEEN  REACHED 

0199 

c 

0200 

IF( DIFMAX. L£ .CTOLIM)  GO  TO  500 

0201 

c 

0202 

c 

INCREMENT  THE  CHORD  COUNTER,  AND  CHECK  TO  SEE 

0203 

c 

WHETHER  THE  MAXIMUM  NUMBER  OF  CHORDS  (MAX  #PT$-1) 

0204 

c 

HAS  BEEN  REACHED. 

0205 

c 

0206 

IF(LINCNT.GT.NOLIN)  GO  TO  500 

0207 

c 

0208 

C' 

NOTE  THE  CHORD  NUMBER  CONTAINING  THE  MAX.  CHORD  DIFF. 

0209 

c 

0210 

LMAX-MAXLIN 

0211 

NMAX-LMAX+1 

j 

j 

i 

i 

i 

t 


?0 


0212 

J=0 

0213 

300 

J=0+1 

0214 

C 

0215 

C 

HAVE  ALL  CHORDS  BEEN  CHECKED 

0216 

C 

0217 

IF(J.GT.LINCNT)  GO  TO  110 

0218 

C 

0219 

C 

DOES  THIS  CHORD  CONTAIN  THE  MAX.  DIFF. 

0220 

C 

0221 

IF(J.EQ.MAXLIN)  GO  TO  120 

0222 

c 

0223 

c 

HAS  THE  CHORD  WITH  THE  MAXIMUM  DIFF.  BEEN  FOUND 

0224 

c 

0225 

lF(J.GT.MAXLIN)  GO  TO  130 

0226 

GO  TO  300 

0227 

120 

K-J+l 

0228 

C 

0229 

C 

ARE  WE  ON  THE  LAST  CHORD 

0230 

C 

0231 

IF(K.EQ.LINCNT)  GO  TO  140 

0232 

c 

0233 

c 

INSERT  THE  NEW  POINT  IN  THE  POINT  PAIR  ARRAY 

0234 

c 

0235 

TEMPl-PTPAlR(K.l) 

0236 

TEMP2«PTPAIR(K,2) 

0237 

PTPAIR(J,2)-MAXPT 

0238 

PTPAIR(K, 1)*MAXPT 

0239 

PTPAIR(K,2)«TEMP1 

0240 

c 

0241 

c 

TEMP.  STORE  THE  DIFFERENCE  FOR  THE  POINT  PAIR  REPLACED 

0242 

c 

IN  THE  ARRAY 

0243 

c 

0244 

DIFTMP-DIFF(K) 

0245 

DIFF(J)«DIFFLT 

0246 

DIFFOO-DIFFLT 

0247 

NUMTMP-NUM(K) 

0248 

NUM(J}*0 

0249 

NUM(K)*0 

0250 

J-J+l 

0251 

GO  TO  300 

0252 

c 

0253 

c 

BUMP  THE  POINT  PAIR,  DIFF.,  AND  MAX.  POINT  ARRAYS 

0254 

c 

DOWN. 

0255 

c 

0256 

130 

TEMP3«PTPAIR(J,1) 

0257 

TEMP4«PTPAIR(J,2) 

0258 

PTPAIR(J ,  1}*>TEMP1 

0259 

PTPAIR(J,2)-TEMP2 

0260 

TEMPI »TEMP 3 

0261 

TEMP2-TEMP4 

0262 

DIFTMP 1«DIFF( J ) 

0263 

DIFF (0 ) *DIFTMP 

0264 

DIFTMP-DIFTMP1 

0265 

NUMTMl-NUM(J) 
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0266 

NUM(J)=NUMTMP 

0267 

NUMTMP=NUMTM1 

0268 

GO  TO  300 

0269 

C 

0270 

C 

TAKE  CARE  OF  THE  LAST  POINTS  IN  THE  ARRAYS 

0271 

C 

0272 

140 

TEMP2-PTPA1R(J ,2) 

0273 

PTPAIR(J,2)-MAXPT 

0274 

PTPAIR(K,1)«MAXPT 

0275 

PTPAIR(K,2)-TEMP2 

0276 

DIFF(J)»DIFFLT 

0277 

DIFF(K)«DIFFLT 

0278 

NUM(J)-0 

0279 

NUM(K)0 

0280 

GO  TO  110 

0281 

500 

NR-LINCNT 

0282 

C 

0283 

C 

USING  THE  POINT  PAIR  ARRAY,  EXTRACT  THE  NEW 

0284 

C 

POINTS  THAT  REPRESENT  THE  THINNED  PROFILE. 

0285 

C 

0286 

L-LINCNT-1 

0287 

DO  2222  0-1, L 

DR(J)-D(PTPAIR(J,1)) 

0288 

0289 

SR(J)-S(PTPAIR(J,1)) 

0290 

2222 

CONTINUE 

0291 

DR(NR)«D(PTPAIR(L,2) ) 

0292 

SR(Nk)-S(PTPAIR(L ,2) ) 

0293 

RETURN 

0294 

END 

0001 

REAL  FUNCTION  INTERP  (P 1 ,Q1 ,P2 ,Q2 , P3 ) 

0002 

C 

0003 

C 

(Pl.Ql)  AND  (P2.Q2)  ARE  PAIRS  BETWEEN  WHICH  OR  PAST  WHICH  THE 

0004 

C 

INTER/EXTRA  -POLATION  IS  DONE.  P3  IS  THE  INTER/EXTRA  -POLATED 

0005 

C 

POSITION. 

0006 

C 

0007 

C 

0008 

IF(P1.NE,P2)  GO  TO  10 

0009 

INTERP «Q1 

0010 

RETURN 

0011 

10 

CONTINUE 

0012 

FACT1-P3-P1 

0013 

FACT2-P2-P1 

0014 

FACT-FACT1/FACT2 

0015 

ANSI -q2-Ql 

0016 

INTERP -F ACT* ANS1+Q1 

0017 

C 

0018 

RETURN 

0019 

END 

0001 

SUBROUTINE  CHORD( YMAX.NUM) 

0002 

PARAMETER  ASIZ=1000 

0003 

COMMON  /DATA/  Y( ASIZ) ,X( ASIZ) 

0004 

COMMON  /PRIM/  XPR1M( ASIZ) ,YPR1M( ASIZ) 

0005 

COMMON  /IPOINT/  1BEG.IEND 

0006 

C 

0007 

c 

DO  AXIS  TRANSLATION  AND  ROTATION 

0008 

c 

0009 

CALL  RTRAN 

0010 

c 

0011 

c 

FIND  THF  MAXIMUM  DEVIATION  OF  THE  TRANSFORMED 

0012 

c 

PROFILE  SEGMENT  FROM  THE  TRANSFORMED  CHORD 

0013 

c 

0014 

YMAX--9999. 

0015 

DO  10  J-IBEG, IEND 

0016 

1F(  ABS(YPRIM(J) ) .LE.YMAX)  GO  TO  10 

0017 

YMAX»ABS(YPRIM( J) ) 

0018 

c 

0019 

c 

NOTE  THE  POINT  NUMBER  OF  THE  MAXIMUM  DEVIATION 

0020 

c 

POINT 

0021 

c 

0022 

NUM-U 

0023 

10 

CONTINUE 

0024 

RETURN 

0025 

END 

0001 

SUBROUTINE  RTRAN 

0002 

C 

0003 

c 

THIS  SUBROUTINE  TRANSLATES  AND  ROTATES  THE  LINE 

0004 

c 

SEGMENTS  SO  THAT  THE  CHORD  BECOMES  THE  NEW  X 

0005 

c 

AXIS.  THE  LAST  POINT  ON  THE  LINE  WILL 

0006 

c 

COINCIDE  WITH  THE  NEW  AXIS  ZERO  (X'-0,Y'-0) 

0007 

c 

0008 

PARAMETER  ASIZ-1000 

0009 

COMMON  /DATA/  Y( ASIZ) ,X(A$IZ) 

0010 

COMMON  /PRIM/  XPRIM{ ASIZ), YPRIM( ASIZ) 

0011 

COMMON  /IPOINT/  IBEG, IEND 

0012 

DO  10  U»IBEG, IEND 

0013 

c 

0014 

c 

TRANSLATE  THE  POINTS  SO  THAT  THE  LAST  POINT  ON  THE 

0015 

c 

LINE  COINCIDES  WITH  0,0  ON  THE  NEW  AXIS 

0016 

c 

0017 

XTEMP*X(J)-X( IEND) 

0018 

YTEMP*Y( IEND)-Y(J) 

0019 

c 

0020 

c 

IF  THIS  IS  THE  FIRST  POINT  OF  THE  LINE  CALCULATE 

0021 

'  c 

THE  ANGLE  OF  ROTATION  REQUIRED  TO  BRING  THE  CHORD 

0022 

c 

DOWN  TO  THE  NEW  X  AXIS. 

0023 

c 

0024 

c 

FOR  THE  FIRST  POINT  XTEMP  IS  DELTA  X  ALONG  THE  LINE 

0025 

c 

AND  YTEMP  IS  DELTA  Y  ALONG  THE  LINE,  WHERE  THE  SLOPI 
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0026 

0027 

0028 

0029 

0030 

0031 

0032 

0033 

0034 

0035 

0036 

0037 

0038 

0039 

0040 

0041 

0042 

0043 

0044 


C  Or  THE  LIME  BETWEEN  THE  TWO  END  POINTS  IS 

C  (DELTA  Y)  /  (DELTA  X). 

C 

IF ( J . EQ . IBEG )  A-ALP HA( XTEMP , YTEMP , 0 . , 0 . ) 

C 

C  FIND  THE  NEW  ROTATED  COORDINATES  FOR  EACH  POINT 

C 

XPRIM(J)-XTEMP*COS( A)+YTEMP*SIN( A) 
YPR1M(U)»-XTEMP*SIN(A)+YTEMP*C0S(A) 

C 

C  IF  THE  LINE  SEGMENT  MUST  BE  ROTATED  BY  90  DEGREES 

C  INSURE  THAT  THE  NEW  Y  COORDINATE  IS  ZERO  AT  THE 

C  END  POINTS. 

C 

IF ( 0 . EO . IB EG. AND. A. EQ. 90.)  YPRIM(0)-0. 

IF( J.EQ. IEND.AND. A.EQ.9Q. )  YPRIM(J)«0. 

10  CONTINUE 
RETURN 
END 


0001 

FUNCTION  ALPHA(X1,Y1,X2,Y2) 

0002 

C 

0003 

C 

THIS  FUNCTION  CALCULATES  THE  ANGLE  OF  ROTATION 

0004 

C 

REQUIRED  TO  ROTATE  THE  CHORD  SO  THAT  IT  LIES 

0005 

C 

ON  THE  X  AXIS. 

0006 

C 

INPUT  ARE  THE  X  AND  Y  COORDINATES  OF  THE  LEFT 

0007 

c 

MOST  POINT  (XI, Yl)  AND  THE  RIGHT  MOST  POINT 

0008 

c 

(X2.Y2). 

0009 

c 

0010 

c 

CHECK  TO  INSURE  AN  INFINITE  RESULT  DOES  NOT  OCCUR 

0011 

c 

0012 

IF(X1.EQ.X2)  GO  TO  1 

0013 

ALPHA-ATAN(  (Y1-Y2)/(X1-X2)) 

0014 

RETURN 

0015 

1 

ALPHA-90. 

0016 

RETURN 

0017 

END 

APPENDIX  C:  VARIABLE  AND  EXTERNAL  CALL  DESCRIPTION 


I.  CHORDS  VARIABLES 

Except  as  noted,  real  variables  begin  with  the  letters  A  through  H,  and  0 
through  Z.  Integer  variables  begin  with  the  letters  I  through  N. 

CTDLIM  -  The  tolerance  value  used  to  halt  processing 

DIFFLT— A  default  value  used  to  fill  unused  portions  of  the  maximum  difference 
array  (DIFF) 

DIFMAX—The  maximum  value  of  all  chord  to  profile  segment  differences 

DIFTMP  and  DIFTMPl--Temporary  storage  locations  for  the  chord  to  profile 
segment  differences  stored  In  array  DIFF 

IBEG--Beg1nn1ng  point  number  of  a  chord  delineated  profile  segment 

IEND--End1ng  point  number  of  a  chord  delineated  profile  segment 

J,K,  and  L--D0  loop  variables  or  Integer  counters 

LINCNT--A  counter  used  to  keep  track  of  the  number  of  chords  found 

LINLIM--A  constant  equal  to  the  total  number  of  chords  allowed 

LMAX--Used  to  store  the  chord  number  containing  the  maximum  chord  to  profile 
segment  difference  of  all  such  sets 

MAXL1N— Same  as  LMAX  above 

MAXPT— The  Input  profile  point  number  defining  the  layer  depth 

MINPT--The  Input  profile  number  point  defining  the  deep  sound  channel  axis 
depth 

NMAX--Used  to  store  the  chord  number  one  greater  than  that  having  the  maximum 
chord  to  profile  segment  difference 

NR— The  number  of  points  In  the  output  thinned  profile 

NUMTM1  and  NUMTMP— Temporary  storage  locations  for  the  number  of  the  chord 
containing  the  maximum  chord  to  profile  segment  difference 

SSMAX—Temporary  storage  location  for  the  value  of  the  Input  profile's  maximum 
above  axis  sound  speed 

SSMIN— Temporary  storage  location  for  the  value  of  the  minimum  Input  profile 
sound  speed 

TEMP3,  TEMP4,  TEMPI,  and  TEMPZ— Integer  temporary  stoage  locations  for  the 
beginning  and  ending  chord  points 


II.  CHORDS  ARRAYS 


D--Input  profile  depth  array 

DIFF  — Used  to  store  the  chord  to  profile  segment  maximum  differences 
DR--0utput  profile  depth  array 

NUM— Used  to  store  the  point  numbers  of  the  maximum  chord  to  profile  segment 
differences 

PTPAIR(N,2)~The  array  containing  the  beginning  chord  point  numbers  In  elements 
(N,l),  and  the  ending  chord  point  numbers  In  elements  (N,2) 

S— The  Input  sound  speed  array 

SR—The  output  sound  speed  array 

III.  CHORDS  REFERENCED  SUBROUTINES  AND  FUNCTIONS 

Function  INTERP(X1,Y1,X2,Y2,X3) 

Interpolates  (or  extrapolates)  between  two  given  points  (XI, Yl)  and 
(X2.Y2)  to  find  the  Y  value  at  X3.  Used  to  determine  the  sound  speed 
at  the  surface  when  It  Is  not  given. 

Subroutine  CHORD (YMAX.NUM) 

This  subroutine  locates  the  maximum  absolute  difference  between  a 
chord  and  Its  corresponding  profile  segment. 


Variables 


IBEG— Beginning  chord  profile  point  number 
IEND--Ending  chord  profile  point  number 

NUM— ' The  number  of  the  profile  point  having  the  maximum  absolute 
difference  from  Its  chord 

YMAX- -Maxi mum  absoulute  difference  between  a  chord  and  Its 
corresponding  profile  segment 


Arrays 


X--Prof11e  segment  depth  values 

XPRIM— Transformed  (translated  and  rotated)  profile  segment  depth 
values 

Y— Profile  segment  sound  speed  values 

YPRIM--Transformed  (translated  and  rotated)  profile  segment  sound 
speed  values 


Subroutines  and  functions  referenced 
RTRAN— (see  below) 

Subroutine  RTRAN 

Transforms  the  sound  speed  profile  segment  points  so  that  the 
corresponding  chord  becomes  the  new  X  axis 

Variables 


A— Angle  of  rotation  of  the  chord  and  profile  segment  such  that  the 
translated  chord  becomes  the  new  horizontal  (X)  axis 

IBEG — Beginning  profile  segment  point 

1END-- Ending  profile  segment  point 

J— Loop  counter 

XTEMP— Temporary  storage  location  of  the  profile  segment's  transformed 
sound  speed  values 

YTEMP— Temporary  storage  location  of  the  profile  segment's  transformed 
depth  values 


Arrays 


X— Profile  segment  sound  speed  values 

XPR1M— Transformed  (translated  and  rotated)  profile  segment  sound 
speed  values 

Y— Profile  segment  depth  values 

YPRIM— Transformed  (translated  and  rotated)  profile  segment  depth 
values 

Subroutines  and  functions  referenced 
Function  ALPHA 

This  function  calculates  the  angle  of  rotation  required  to  rotate  the 
chord  so  that  It  becomes  parallel  to  the  X  axis 
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